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ABSTRACT 

We have considered linear two point correlations of the form 4^ which are known to 
have a self-similar behaviour in a 17 = 1 universe. We investigate under what conditions 
the non-linear corrections, calculated using the Zel'dovich approximation, have the 
same self-similar behaviour. We find that the scaling properties of the non-linear 
corrections are decided by the spatial behaviour of the linear pair velocity dispersion 
and it is only for the cases where this quantity keeps on increasing as a power law 
(i.e. for 7 < 2) do the non-linear corrections have the same self-similar behaviour as 
the linear correlations. For (7 > 2) we find that the pair velocity dispersion reaches 
a constant value and the self-similarity is broken by the non-linear corrections. We 
find that the scaling properties calculated using the Zel'dovich approximation are very 
similar to those obtained at the lowest order of non-linearity in gravitational dynamics 
and we propose that the scaling properties of the non-linear corrections in perturbative 
gravitational dynamics also are decided by the spatial behaviour of the linear pair 
velocity dispersion. 



Subject headings: Galaxies: Clustering - Large Scale Structure of the Universe 
methods: analytical 



1. Introduction. 



The equations governing the evolution of the statistical properties of disturbances in a critical 
(f2 = 1), matter dominated universe are known to admit self-similar solutions (Peebles 1980). 
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This is because the universe expands as a power law of time i.e. a(t) oc ts and gravity itself does 
not introduce any preferred scale. As a consequence of this it is possible, using the linear theory 
of density perturbations, to construct correlation functions which have a self-similar behaviour 
over a range of scales. This method is valid only as long as the correlations on these scales are 
extremely small. On the other hand it is possible to use the stable clustering assumption to 
construct correlation functions that have a self-similar behaviour over a range of scales. The 
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latter assumption is valid only on small scales where virialized objects have already formed and 
clustering in real space has ceased to increase. On scales where neither of these assumptions can 
be used very little is known. 

Hamilton et. al. (1991) have suggested a universal scaling relation for the two point 
correlation function based on numerical evidence from N-body simulations. In more recent papers 
Nityananda & Padmanabhan (1994), Peacock & Dodds (1994), Jain et. al. (1995), Lokas et.al. 
(1995) and Padmanabhan (1996) have investigated the proposed universal scaling relations. In 
an earlier paper (Bharadwaj 1996) we have studied the lowest order non-linear corrections to the 
two point correlation function for cases where the linear correlation function has a self-similar 
behaviour over a range of scales. We restricted ourselves to cases where the linear power spectrum 
has the form P(k) oc k n with n > at small k and we found that the non-linear corrections do not 
have the same self-similar behaviour as the linear two point correlation function. In a recent paper 
Scoccimarro and Frieman (1996) have addressed the same question for initial conditions which 
include n < 0. For n > their conclusions are similar to ours, but they find that for n < —1 
the lowest order non-linear correction has the same self-similar behaviour as the linear power 
spectrum. 

In this paper we consider initial perturbations which are a random Gaussian field and have a 
two point correlation function which is self-similar over a range of scales. We study the evolution 
of the correlation function analytically using the Zel'dovich approximation (ZA) (Zel'dovich 1970) 
and we investigate under what general conditions the non-linear corrections to the two point 
correlation function exhibit the same self-similar behaviour as the linear two point correlation 
function. ZA is known to be a very good approximation to the full gravitational dynamics in 
the weakly non- linear regime before the effects of multi-streaming become important. It is thus 
expected that the results of this investigation should also hold for gravitational dynamics and we 
have compared our results with the results available from perturbative non-linear gravitational 
dynamics (GD). The use of ZA instead of GD makes the algebra much more tractable and it is 
much simpler to interpret the results. It is also hoped that the results of this investigation may 
help in building up correlation functions that are self-similar in the non-linear epoch starting from 
correlations that are self-similar in the linear regime. 

The evolution of the various correlation functions in ZA has been studied by earlier authors 
including Bond & Couchman (1988), Grinstein & Wise (1987), and Schneider Sz Bartlemann(1995 
). In this paper we use the notation and some of the results of our earlier paper (Bharadwaj 1995) 
where we also presented some of our preliminary results regarding the non-linear corrections to 
the two point correlation function in ZA and a comparison with the results from GD. 



2. The two point correlation in ZA. 
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For a ft = 1 universe ZA defines a map 

x tl (t) = x fl (to) + a(t)u tl (1) 

from the initial position x M (io) of a particle to its position x M (i) at some later time t in a comoving 
coordinate system. Here the subscript /i takes values 1,2 and 3 corresponding to the three 
Cartesian components and the Einstein summation convention holds for it. The quantity u M is 
related to the peculiar velocity of the particle and we shall refer to it as the velocity. 

We are interested in the evolution of the statistical properties of an ensemble of systems whose 
evolution is governed by ZA. In all the members of the ensemble the particles are all initially 
uniformly distributed and the initial velocity field is assumed to be irrotational. It is also assumed 
that the velocity field in any member of the ensemble is a particular realization of a Gaussian 
random field. These initial conditions can be fully specified by the velocity-velocity correlation 
which can be written in terms of a potential <f>{x) as 

< u^x^u^x 2 ) >= -d^d u (f>(x) (2) 



where dp = and 



x =\ x — x | . (3) 



For such an ensemble the two point correlation function can be written in a perturbative 
expansion as (Bharadwaj 1995) 

1 



d^d^x) - ^ m V=V(0) 



£(M) = Y,™ = y n -d^d vl d^d V2 ...d^d Vn 

\ 2 d„ 2 cj>(x) - ^<W 2 V 2 <K0)) ... [d^d^x) - I^ n V 2 <K0) 



(4) 

The quantity that appears in the right hand side of this equation is the dispersion of the pair 
velocity (i.e. the difference between the velocity at the point x 2 and the point x l ) which is defined 

as 

< v^Vy > (x) = < (u^(x 2 ) - u M (x 1 ))(u M (x 2 ) - ^(x 1 )) > 

= 2 (d^d^ix) - ^<W 2 V 2 <M0)) . (5) 

This quantity which we shall refer to the as the pair velocity dispersion is related to the dispersion 
of the relative peculiar velocities a^ v (x) in the linear epoch, and we have 

da(t)\ 2 



cr^u(x,t) = (^W - ^ - J < v v v » > ( x ) ( 6 ) 

This dispersion arises due to the spread in the relative velocities across the various realizations 
in the ensemble. The pair velocity dispersion is a symmetric tensor and because the initial 
conditions are statistically homogeneous and isotropic, it can in general be written as 

< v„v v > (x) >= 5, v P(x) + ^Q(x) (7) 
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where P{x) is the dispersion of the relative velocity component perpendicular to the separation 
x(= x 2 — x 1 ) and P(x) + Q(x) is the dispersion in the velocity component parallel to x. 

The two point correlation function can be written in terms of the pair velocity dispersion as 

^ ( a 2 \ n 

C(x,t) = ££° =1 — f y J d^d^d^d^.-.d^d^ [< v^v Vl > (x) < v^ 2 v V2 > (x) 

- < v„ n v Un > (x)} . (8) 

From this we obtain the linear two point correlation function as 

C (1 \x,t) = a 2 ^(x) = y^Ax < Vlxl v vl > (x) = a 2 VV(x) . (9) 

Given the linear two point correlation function we can invert equation (|9|) and express the 
initial pair velocity dispersion in terms of the linear two point correlation function as 

< Vfj,v v > (x) = -8^ / i w {y)ydy 

O JO 

x 







-d^x) I i {1) (y)y 2 dy 

1 

3 



d»d v {l/x) / ^ 1) {y)y A dy (10) 
J o 

and finally we are in a position to start investigating the nature of the non-linear corrections to 
the two point correlation function, given the linear two point correlation function. 

We consider situations where the linear two point correlation function at some large scales i.e. 
for x > x\ has the form 

{«(«,,) _•>(!) (ay _(2WV (ii) 



where 

2 / t \ ^ 

xo(t) = x a(t)i = x ^— J . (12) 

We see that the two point correlation function has a self-similar behaviour over a range of scales 
in the linear epoch and the effect of temporal evolution is to just scale the length scale xo(t) that 
appears in the two point correlation function. We want to investigate under what conditions the 
non-linear corrections to the two point correlation function has the same self-similar behaviour. 
Since the non-linear corrections are determined by the pair velocity dispersion, we proceed by first 
investigating the behaviour of the pair velocity dispersion. 



3. The pair velocity dispersion. 
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In this section we investigate the spatial behaviour of the pair velocity dispersion for different 
values of the index 7. We separately consider the three different moments of the linear two point 
correlation function that appear in equation (pX|) for the dispersion of the pair velocity 

A(x)= ! X ^\y)ydy (13) 
Jo 

B(x) = - r^Ky)y 2 dy (14) 
x Jo 

C(x) = ± f x ^ 1 Hy)y i dy (15) 

x 6 Jo 

It should be noted that all three of the functions defined above have dimension I? and it is 
possible that any of them may introduce a new length scale in the evolution. 



3.1. The behaviour of A. 



For 7 > 2 the integral in A converges in the limit x — > 00 and at large x we have 

A(x) = / e ) (y)ydy+( _2_) (16) 







2-7/ XT- 2 • 



Here the first term is a constant and the second term decays as x increases, and at large x A(x) 
tends to a constant value.. 

At the value 7 = 2 we have 

A[x) = jH (y)ydy + xl In (£) (17) 



For 7 < 2 we have 



where the first two terms are constants and the last term increases monotonically with x and for 
large x we have 

A & = (2^) ^ (19) 



3.2. The behaviour of B 
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For 7 > 3 the integral in B converges as x — > oo and at large x we have 



1 roo 

B(x) = - t {l) (y)y 2 dy 
x Jo 



+ 



3 - 7 J 



1 



(20) 



Here the first term falls as ^ and the second term falls off faster than this, and hence at large x 
the first term dominates. 



At the value 7 = 3 we have 



B{x) 



and for large x we have 



/ \ {1) {y)ydy + 
Jo 

B(x) 



3 

Xq fx 
— log 

X 



and for 7 < 3 we have 

B(x) = -( re\y)y 2 dy 
x \Jo 



3-7 



x 



3-7 



+ 



3-7 



x 



2-7 



(21) 
(22) 

(23) 



where for large x the first terms falls off as ^ and the second term falls off slower. Thus at large x 
we have 



B{x) = 



3-7 



9_ \ 



X 



(24) 



where this is a decaying function for 2 < 7 < 3, it is a constant for 7 = 2 and it is a monotonically 
increasing function for 7 < 2. 



3.3. The behaviour of C 



For 7 > 5 the integral in C converges asx^oo and at large x we have 



1 r°° 

c(x) = - ^(y)y 4 dy 



+ 



5 — 7 / x~< 



(25) 



here the first term falls as ^ and the second term falls off faster than this, and hence at large x 
the first term dominates. 

At the value 7 = 5 we have 

a,) = J 

and for large x we have 



Xl t {l) {y)y 4 dy + xl\n (- 
\xi 



C(x) 



(26) 
(27) 
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For 7 < 5 we have 

= ^ (f - 4-) + (^) (28) 

where at large x the first terms falls off as and the second term dominates. Thus at large x we 
have / 7 \ 

C(x) = (^J . (29) 

where this is a decaying function for 2 < 7 < 3, it is a constant for 7 = 2 and it is a monotonically 
increasing function for 7 < 2. 



3.4. The spatial behaviour of the dispersion of pair velocities. 

Putting together the various components calculated earlier we can now write expressions for 
the dispersion of the pair velocity for different values of the index 7. 

For the cases where 7 < 2 the behaviour is very simple and we have 

< v,v v > (x) = V (^^y " (3^) + ^^y) *> 2 ~ 7 

^ ( — —) x2x 2 ^ 

x \ 3 — 7 5 — 7 / 



+ ( U> 7 (30) 



For this case we see that both the components of the pair velocity dispersion increase as oc x 1 7 . 
Next, for 7 = 2 we get 



< VpVp > (x) = ^Xq 



(31) 



Here both the components of the pair velocity dispersion increase as ln(x) and are nearly equal as 
we go to large separations 

For 2 < 7 < 3 we have to add an extra term 



•J /'IAJ 

y 2 = ^( t (l \y)ydy (32) 







to equation fl30|). This term dominates the behaviour of the pair velocity dispersion at large 
x as the contribution from the terms in equation ( |30| ) get smaller as x increases and both the 
components of the pair velocity dispersion tend to a constant value l 2 . Thus we see that we have 
one new length scale now i.e. I and this plays a very crucial in the later discussion. 
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For 7 = 3 we have 



< VpVv > (x) 



9 roo ~3 

-J»v e (1) (y)y^-V- 

3 Jo x 



1 , fx 
- + In — 

2 \xi 



. t p dj V J Q 



(33) 



For 3 < 7 < 5, in addition to the expression (p2|), we have to also add 

^ l \y)y 2 dy 



Xy,X v j- 








.a; 


x Jo 



(34) 



to equation (30). 

For 7 = 5 we have 



2 Z" 00 

>W = A £ (1) (y)y*/ + 

o Jo 



X p^X v 



J fj,u 



+5 -^2. 

+V 3x3 



5 , / a; 
- + In — 

6 V^l 



x Jo 

,5 



^\y)y 2 dy 



X 2 X 3 



For 7 > 5 , in addition to the two terms in expressions ( |32|) and 

-/ £ (1) (y)y 4 */ 

X 13 Jo 



1 , fx 
- + In — 

2 Ui 



we have to also add 



1 r XyJX v 

3 V " -^2" 



(35) 



(36) 



to equation (30) 



Finally, we see that for 7 < 2 the dispersion in the pair velocity keeps on increasing as a power 
law. For 7 = 2 it increases logarithmically with the separation and for all 7 < 2 pair velocity 
dispersion tends to infinity as the separation x tends to infinity. For 7 > 2 we find that the pair 
velocity dispersion reaches the constant value I 2 which is determined by the two point correlation 
at small separations. 

We next consider some specific cases where the linear power spectrum is of the form 
P(k) = Ae~ k k n with n = (a) — 2, (b) — 1, (c) and (d) 1. The corresponding values of 7 are 
(a) 1, (b) 2, (c) 4 and (d) 4. For the case with n = the large x behaviour is decided by the 
exponential cutoff in P(k) and because of this it has the same power law index (7 = 4) as the 
n = 1 case. Figure (1) shows the tangential component of the pair velocity dispersion for these 
four cases and it illustrates the main point of this section. The behaviour of the radial component 
of the pair velocity dispersion is very similar. 



4. The non-linear corrections to £. 
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Here we use the pair velocity dispersion calculated in the previous section to investigate the 
nature of the non-linear corrections to the two point correlation function. From equation (|l^) we 
obtain the lowest order non-linear correction to the two point correlation function as 



— d^d^d^d^ [< v 



^iVui > (x) < v^ 2 v V2 > (x)] . 



(37) 



For 7 < 2 the lowest order non-linear correction to the two point correlation function is 
2(285 - 6797 + 6II7 2 - 2597 3 + 52 7 4 - 4 7 5N i /t«\ 27 



e (2) (; 



(2- 7 )(3-7) 2 (5-7) S 



a-oY 

X J 



(38) 



and for this range of 7 we find that £ (2) (M) oc £ (1) (a;,t) • For 7 < 2 the pair velocity dispersion 
is a power law in x and it is ~ x 2 ~ 7 . To calculate the n th order non-linear corrections to the two 
point correlation function we take n of these functions i.e. an expression cx (x 2 ~ 1 ) n and we act on 
this with 2n spatial derivatives and the result is x~ ni oc >{x) . Thus we see that nth order 
non-linear correction is of the form ^ n \x) oc [^^(x)]™ and we can write the two point correlation 
function in a series of the form 



(39) 



where the c n (^) are coefficients which depend only on 7. We see that for 7 < 2 the non-linear 
corrections to the two point correlation function have the same self-similar behaviour as the linear 
two point correlation function. 



(40) 



For 7 = 2 we have 

For this case we find that the non-linear correction cannot be written in terms of the linear two 
point correlation alone and it does not have a self-similar behaviour. 



These are the only cases where the pair velocity dispersion keeps on increasing at large x. 
For larger values of 7 the pair velocity dispersion reaches a constant value and as consequence the 
behaviour of the two point correlation too is different. For all higher values of 7, in addition to 
the contribution to C^ 2 \x) given by equation ( |38| ) which behaves as x^ 4_27 \ there are other terms 
which fall of slower and dominate the behaviour at large x. For (2 < 7 < 3) we have 



(2), 



^{y)ydy 



7(7 - 1)- 



(41) 



For these cases we see that the self-similar behaviour is broken. This is because the pair velocity 
dispersion reaches a constant value and this introduces a new length scale I whose square appears 
in the square brackets in equation (|4l]) . This length scale I changes the scaling property of the two 
point correlation function. 
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For 7 > 3 we have to take into account one more term and we have 



£ (2) M) 



1 r°° 1 t" 1 3 r r°° 

- e\y)ydy 7 ( 7 _l)^- + _ / e\y)y 2 dy 

o JQ X ! X \_J0 



2 



(42) 



For 7 < 4 it si obvious that the first term falls off slower than x~ 6 and it dominates the large 
x behaviour. For 7 = 4 both the terms have a x -6 behaviour and we cannot drop any one of 
them. For 7 > 4 the first term falls off faster than x~ 6 and we would expect the second term to 
dominate, but for many of the cases of interest the integral f£° i <yl \y)y 2 dy is zero and then it is 
only the first term that contributes. This is illustrated by the two cases (c) and (d) that we have 
considered earlier. The moments of the two point correlation function can be written in terms of 
the power spectrum as 

t (1) (y)ydy = ^ J q P(k)dk (43) 

and 

'e\y)y 2 dy = ^. (44) 

47T 

For both (c) and (d) we have 7=4, but using the above equations we find that for (c) (i.e.n = 0) 
we obtain J °° ^ 1 \y)y 2 dy = ■S- and for (d) (i.e.n = 1) we have / °° i^ l \y)y 2 dy = 0, and for both of 
them we have Jq°° ^ l \y)ydy = j^. Thus for (c) the second term in equation (f42|) contributes and 
for (d) it is zero and for these two cases we obtain 

to (a) ^ X ) = ^( 2x ° + ^)^ («) 

and 

for(b) ^ 2) (x) = ^l. (46) 

7T Z X D 

Finally we can generalize this to say that for all cases where linear power spectrum is of the 
form P(k) oc k n with n > — 1 for small k, the lowest order non-linear correction to the two point 



correlation function is given by equation (41), except for the case when n = 0. For n = we have 



to take into account the extra term in equation ([421). 

In a recent paper Taylor and Hamilton (1996) have considered the evolution of the non-linear 
power spectrum in the Zel'dovich approximation. They present exact analytic results for cases 
where the linear power spectrum is a power law with the power law index n = —2, —1 i.e. (7 = 1, 2). 
We have compared our result at the lowest order of non-linearity with the corresponding result 
implied by the exact expression calculated by Taylor and Hamilton (1996). We find that while the 
two match for n = —2, there is a disagreement for n = —1. This is because for n = —1 (7 = 2) 
they have considered the pair velocity dispersion to bee a constant value whereas we find that it 
has a logarithmic dependence on the separation. In addition to this, they have considered n = —1 
as a limiting case of a situation where the power law index is of the form n = — 1 — e. In our study 
here we find that the behaviour of the case with 7 = 2 i.e. n = — 1 is quite different from the 
results for 7 < 2 and we do not expect that the limit taken by Taylor and Hamilton (1996) will 
give the correct result for n = — 1. 
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5. Discussion and Conclusion. 

We have considered linear two point correlation functions that have the form ^ at large 
separations. For = 1 they have a self-similar behaviour in the linear epoch. We have investigated 
under what conditions the non-linear corrections calculated using the Zel'dovich approximation 
have the same self-similar behaviour. 

We find that the scaling properties of the non-linear corrections to the two point correlation 
function are determined by the spatial behaviour of the linear pair velocity dispersion. For 7 < 2 
both the radial and tangential components of the pair velocity dispersion keep on increasing 
as a; 2-7 and it has a local dependence on the linear two point correlation function i.e. the 
linear pair velocity dispersion at some separation x depends only only the linear two point 
correlation at the same separation. As a consequence all the non-linear corrections to the two 
point correlation function also have a local dependence on the linear two point correlation function 
and the n th order non-linear correction to the correlation function £} n \x,t) has the property 
^ n \x,t) ~ [^{x,t)] n . We see that for these cases all the non-linear corrections have the same 
self-similar behaviour as linear two point correlation. 

For 7 = 2 the pair velocity dispersion increases logarithmically at large x and the n th order 
non-linear correction is of the form ^^%- f — i.e. ^ n \x,t) ~ [In ( — )]" -1 [^ 1 Hx, t)] n . Because of 

x y x x J 

the extra [In (j^] n ~ l the non-linear corrections do not have the same self-similar behaviour as 
the linear two point correlation function. 

For 7 > 2 the pair velocity dispersion at large x reaches a constant value I 2 which is 
determined by the linear two point correlation at small scales. For this case the linear pair 
velocity dispersion has a non-local dependence on the linear two point correlation function. This 
also introduces a new length scale and the n th order correction is of the form * + ( 2 (n-i) ■ This 
has completely different scaling properties and the non-linear corrections do not have the same 
self-similar behaviour as the linear two pint correlation. 

In an earlier paper (Bharadwaj 1996) we have studied the non-linear correction to the 
two point correlation function using perturbative gravitational dynamics. We only considered 
cases with 7 > 2 and we found that the non-linear corrections do not have the same self-similar 
behaviour as the linear two point correlation function. We also found that this was due to the 
emergence of a new length scale I which appears here also and we interpreted this in terms of a 
simple diffusion process. Scoccimarro & Frieman (1996) have studied the lowest order non-linear 
correction to the power spectrum and the corresponding correlation functions include those with 
7 < 2. They find that for power spectra where the index n is less than —1 (or 7 < 2), the lowest 
order non-linear correction has the same self-similar behaviour as the linear power spectrum and 
they also find that the self-similarity is broken for n > — 1 (i.e. 7 > 2). We see that ZA makes 
the same predictions as GD regarding the scaling properties of the non-linear corrections to the 
two point correlation function and we expect that any conclusions that can be drawn on the 
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basis of our investigations using ZA should also hold for GD. We find that in ZA the scaling 
properties of the non-linear corrections are decided by the spatial behaviour of the linear pair 
velocity dispersion. We propose that this is also true for the non-linear corrections calculated using 
perturbative gravitational dynamics. Although the relation between the non-linear corrections 
and the linear pair velocity dispersion is quite explicit in the Zel'dovich approximation (equation 
H), it is not clear how this comes about in perturbative gravitational dynamics. 

The pair velocity dispersion is the lowest velocity moment of the two point distribution 
function that has information about the velocity components tangential to the separation x. It is 
very interesting that it is this quantity, and not some lower velocity moment, that turns out to 
play an important role in deciding the scaling properties of the non-linear corrections. 

The author would like to thank Rajaram Nityananda for many useful discussions and T.R. 
Seshadri for his encouragement. 
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Fig. 1. — The tangential component of the pair velocity dispersion P(x) shown as a function of the 
separation x for the four cases discussed in section 3.4. 



